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In this paper, we investigate an optimal investment and con- 
sumption problem for an investor who trades in a Black-Scholes fi- 
nancial market with stochastic coefficients driven by a non-Gaussian 
Ornstein-Uhlenbeck process. We assume that an agent makes invest- 
ment and consumption decisions based on a power utility function. 
By applying the usual separation method in the variables, we are 
faced with the problem of solving a nonlinear (semilinear) first-order 
partial integro-differential equation. A candidate solution is derived 
via the Feynman-Kac representation. By using the properties of an 
operator defined in a suitable function space, we prove uniqueness 
and smoothness of the solution. Optimality is verified by applying a 
classical verification theorem. 

1. Introduction. A fundamental problem in financial mathematics is the 
allocation of funds between assets in order to provide sufficiently large pay- 
ments during the duration of an investment contract, as well as to arrive 
at a high return at maturity. This optimization problem has its origin in a 
seminal paper by Merton [18], where it is formulated as a utility maximiza- 
tion problem and an optimal strategy is derived via the Bellman equation. 
Since then, there has been a growing interest in investment and consumption 
problems and the classical Merton problem has been extended in many di- 
rections. One of the generalizations considers financial coefficients (risk- free 
return, drift and volatility) affected by an external stochastic factor. 
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In this paper, we extend the results from [5] and [17]. We investigate a 
Black-Scholes-type financial model with coefficients depending on a back- 
ground driving process. The dependence is described through general func- 
tions which satisfy linear growth conditions. An external stochastic factor 
is chosen as Ornstein-Uhlenbeck process driven by a subordinator. The 
Barndorff-Nielsen and Shephard model considered in [5] and [17] arises as a 
special case. As an additional possibility, the investor is allowed to withdraw 
(consume) funds during the term of the contract. This leads to an optimal 
investment and consumption problem which is more complex than a pure 
investment problem. From an analytical point of view, the difference is that, 
after applying the usual separation of variables, we arrive at a nonlinear 
partial integro-differential equation, whereas [5] and [17] deal with a linear 
one. 

The first goal of this paper is to show that a candidate value function 
is the classical solution of a corresponding Hamilton-Jacobi-Bellman equa- 
tion. This requires proving the existence of a classical solution to a nonlinear 
(semilinear) first-order partial integro-differential equation. It is well known 
(see [9], Chapter 12.2, and [20]) that the regularity of solutions to equa- 
tions with an integral term is uncertain, especially in the degenerate case. 
There exist some results concerning the smoothness of a solution to a lin- 
ear partial integro-differential equation (see, e.g., [2], Chapters 3.3 and 3.8, 
[9], Chapter 12.2 and [20]), but they all deal only with the nondegenerate 
second-order case. The degenerate case can be handled by applying a vis- 
cosity approach (e.g., [9], Chapter 12.2) which we want to avoid, following 
instead [5], where the existence of a classical solution to a linear first-order 
partial integro-differential equation is established. We believe that our proof 
(in Sections 4 and 5) of the existence of a unique classical solution to a 
nonlinear first-order partial integro-differential equation contributes to the 
present state of the literature. 

Our second goal is to provide an explicit formula for the optimal consump- 
tion. In the case of a power utility function, it is intuitively easy to foresee a 
formula for the optimal investment, by simply replacing deterministic coef- 
ficients by functions, which relate coefficients to an external factor and thus 
adapt the strategy to an underlying filtration. This is no longer obvious as 
far as the consumption strategy is concerned. To the best of our knowledge, 
the formula for the optimal consumption in the model investigated in this 
paper is new (see Theorem 6.1). 

Portfolio optimization in stochastic factor models has recently gained 
much attention in the financial literature. In the majority of papers, a power 
utility function is applied and a Black-Scholes financial market with an ex- 
ternal stochastic factor of diffusion type is considered. In this setting, it is 
well known that one must solve a nondegenerate nonlinear second-order par- 
tial differential equation. Several methods have been proposed to deal with 



OPTIMAL INVESTMENT AND CONSUMPTION 



3 



this problem. In [23], in the case of a pure investment problem, a power 
transformation was introduced, which makes the nonlinear term disappear. 
In [16] a similar transformation has been applied, but because of the pos- 
sible consumption, a linear partial differential equation appears only in the 
case of perfectly (positively) correlated Brownian motions or for logarithmic 
utility. More effective methods have been proposed in [12] and [8]. In the 
first paper, a change of measure transformation is applied and the resulting 
optimization control problem is investigated, whose value function depends 
only on time and a factor variable. In the second paper, the dual problem is 
considered, whose control process belongs to a set of equivalent local mar- 
tingale measures. Again, the value function of the dual problem depends 
only on time and a factor variable. This method has also been successfully 
applied in a robust utility maximization model in [15] recently. In all three 
aforementioned papers, the existence of a classical solution to the Hamilton- 
Jacobi-Bellman equation is proved in three steps: first, by constraining the 
values of the control process to a compact set, second, by applying results 
from the theory of nondegenerate linear partial differential equations (see 
Chapter VI. 6 and Appendix E in [13]), thus showing that the constrained 
problem has a unique classical solution, and, third, by studying the asymp- 
totic limit. It seems that this method cannot be successfully applied to our 
problem. 

In the present paper, a candidate solution is first derived heuristically via 
the Feynman-Kac representation. This leads to a fixed point equation. The 
existence of a solution is established by Banach's fixed point theorem and 
its differentiability is proved by using the properties of a suitable operator. 
Finally, we show that the candidate solution satisfies our integro-differential 
equation and that this solution is unique. The idea of finding a solution to a 
control problem through a fixed point theorem is not new; it is, for example, 
mentioned in [8]. In [4], the existence of a solution to a nondegenerate non- 
linear (semilinear) partial differential equation is proved by Banach's fixed 
point theorem. The smoothness then follows from Holder estimates for a 
solution of a nondegenerate linear partial differential equation. We would 
like to point out that, in particular, in [4], an exponent in the Feynman-Kac 
formula is assumed to be bounded, which leads to a bounded solution, while 
we are dealing with a solution which satisfies only an exponential growth 
condition. We would also like to mention that in the context of optimal con- 
trol, the results from [4] are directly applied in [10], where an investment and 
consumption problem is investigated in the presence of default, triggered by 
a one-jump counting process with a stochastic intensity of diffusion type. 

Throughout this paper, we assume that the external factor is observable 
(as in all aforementioned publications). An alternative would be a partially 
observed control problem, whose optimal strategy would then be based on 
an estimate of the underlying factor. We refer to [3] or [21], where a portfolio 
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problem is solved in a diffusion setting with an unobserved volatility process 
of diffusion type and of Markov switching type, respectively. 

Our paper is structured as follows. In Section 2, we introduce the fi- 
nancial market. The optimization problem is formulated in Section 3. The 
uniqueness of a solution is proved in Section 4, whereas the differentiability 
is established in Section 5. In Section 6, we show the optimality of a solu- 
tion and illustrate our findings by means of a numerical example. We also 
present the solution to the optimal investment and consumption problem 
for logarithmic utility. 

2. The financial market. Let (f2, J 7 , IP) be a probability space with filtra- 
tion F= (J-(t))o<t<T, where T denotes a finite time horizon. The filtration 
is assumed to satisfy the usual conditions of completeness and right conti- 
nuity. The measure P is the real- world, objective probability measure. All 
expectations are taken with respect to P. 

We consider a Black-Scholes market with coefficients driven by an ex- 
ternal stochastic factor. Let Y := (Y(t))o<t<T denote this economic factor, 
whose dynamics is given by a stochastic differential (SDE) equation of the 
Ornstein-Uhlenbeck type, 

(2.1) dY(t) = -\Y(t-)dt + dL(Xt), Y(0) = y>0, 

where A > denotes the reversion rate and L := (L(t))o<t<T is an F-adapted 
subordinator with cadlag sample paths. Recall that a subordinator is a Levy 
process with a.s. nondecreasing sample paths. For definitions and more back- 
ground on Levy processes, we refer to [1, 7] or [22]. 

Our financial market consists of two instruments. The price of a (locally) 
risk- free asset B := {B(t))o<t<T is described by the differential equation 

/ x dB(t) , , „ 

(2.2) ^L =r (Y(t-))dt, 5(0) = 1, 

whereas the dynamics of the price of a risky asset, S := (S(t))o< t <T, is given 
by the SDE 

(2-3) ^ = n{Y (t-)) dt + a(Y(t-)) dW(t), 5(0) = s > 0, 

where W := (W(t))o< t <T denotes an F-adapted Brownian motion, indepen- 
dent of the subordinator L. We make the following assumptions concerning 
the functions r,fx and a: 

(Al) the functions r:(0, oo) — > [0, oo), /x:(0, oo) — ► [0, oo) and a : (0, oo) — > 
(0, oo) are continuous and satisfy the linear growth conditions 

r(y) <A r + B r y, fi(y) < + B^y, 

a 2 (y)<A a + B a y, y>0, 

with nonnegative constants; 
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(A2) the derivatives J : (0, oo) -> K, ^ : (0, oo) -> K and ^ : (0, oo) -> R are 
continuous and satisfy linear growth conditions analogous to those of 
r,/x,cr 2 ; 

(A3) infygDa f(y) > 0, where the set T>2 will be specified in (3.7). 

Note that the assumptions (A1)-(A3) are more general than in [8, 12] 
and [15], where uniform boundedness of the functions r, /i, a 2 and their first 
derivatives is required. Our conditions are similar to those in [23], where 
Lipschitz continuity of the coefficients is assumed, together with a linear 
growth condition. 

A prominent example of the above financial model is the Barndorff-Nielsen 
and Shephard model, introduced in [6], which can be described by the fol- 
lowing set of equations: 

dB{t) 



Bit) 
dS(t) 



rdt, 



(2.4) — ^ = (,i + 0Y(t-)) dt + yjY(t-) dW(t). 

Besides the above paper, we also refer to [9], Chapter 15, [5, 17] and ref- 
erences therein for more information about the properties of non-Gaussian 
stochastic volatility models in the context relevant to our paper. 

We shall need some further results and notation for Y and its background 
driving Levy process L. The subordinator L has the representation (see, e.g., 
[1], Chapter 1.3.2) 



(2.5) L(t)= / / zN(ds,dz), t > 0, 

J(P,t] Jz>0 

where N((0,t] x A) = #{0 < s <t:(L(s) - L(s-)) G A} denotes a Poisson 
random measure with a deterministic, time-homogeneous intensity measure 
v(dz) ds satisfying / <^<i Zl/ (dz) < oo. The fundamental result in the theory 
of infinitely divisible random variables is the Levy-Kintchine formula, which 
presents the moment generating function of a subordinator as 



w <w, 



(2.6) E[e wL ^}=e t ^ w) =exp^tj^ (e wz - l)v(dz)\, 

for some w E [0, oo]. The function ipiw) is called the Laplace exponent of L. 
Note that tp(w) exists at least for all w < and tp(w) > for w > 0, provided 
it exists. 

Let us now investigate the SDE (2.1). Its unique solution for s > t is given 
by (cf. [1], Chapter 6.3) 

(2.7) Y(s) = ye- x{s - t) + f S e-^-^dLiXu), Y(t) = y. 



(2.8) -K-Y t,v ( 



-A(s-t) 



s = e v ' r-a.s. 
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We abbreviate the process (2.7) by Y t,v := (Y t ' y (s),t <s<T) and would 
like to point out that it has a.s. cadlag sample paths of finite variation and 
that the mapping y i— > Y t,y is continuous P-a.s. Moreover, 

d_ 

dy 

Finally, it is straightforward to establish the following relations for all < 
t < s < T and y > 0: 

(2.9) Y^(s)<y + L(Xs)-L(Xt), 

(2.10) A J°Y t < y {u)du = y + L(Xs)-L(Xt)-Y t < y (s) 

(2.11) < y + L(Xs) - L(At) = y + L(X(s - t)). 

The above relations hold P-a.s., except for the last equality, which holds in 
distribution. 

3. Formulation of the optimization problem. We consider an investor 
who makes decisions concerning investment and consumption of a portfolio 
based on a power utility function of the form for 7 G (0, 1). 

Consider the wealth process X C,1T := (X c,7r (t))o<t<T of an agent. Its dy- 
namics is given by the stochastic differential equation 

dX c ^(t) = Tr(t)X c ^(t)(n(Y(t-)) dt + a(Y(t-)) dW(t)) 

(3.1) 

+ (1 - TT{t))X c ' n {t)r(Y{t-)) dt - c{t) dt, 

where ir(t) denotes a fraction of the wealth invested in the risky asset and 
c(t) denotes the rate of consumption at time t. We are dealing with the 
following optimization problem: 



(3.2) supE 



T 







(c(s)r ds + (X^(T)y I X(0) = x, Y(0) = y 



The corresponding optimal value function is defined as 
V(t, x, y) = sup E 

(3.3) 



(c( s )yds + (x c ^(T)y\x(t) = x, 

Y(t) = y 



Let us introduce the set A of admissible strategies. 

Definition 3.1. A strategy (c,tt) := (c(i),7r(i))o<t<r is admissible, and 
we write (c, ir) £ A, if it satisfies the following conditions: 
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1. (c,7r) : (0,T] x — > [0, oo) x [0, 1] is a progressively measurable mapping 
with respect to the filtration F; 

2. Jq c(s) ds < oo P-a.s.; 

3. the SDE (3.1) has a unique, positive solution X C,7T on [0,T]. 

We would like to mention that for every (c, ir) € A, the wealth process 
X c ' n , which satisfies (3.1), is an Ito diffusion; that is, in particular, a semi- 
martingale with P-a.s. continuous sample paths. 

Note that we exclude the possibility of borrowing from the bank account 
and short-selling the asset, as in [5] and [17]. Technically, there is no problem 
in solving the unconstrained optimization problem. In particular, if (fi(y) — 
r(y))/(T 2 (y) is positive and uniformly bounded, then all of our results remain 
the same. 

One can associate a Hamilton-Jacobi-Bellman equation with the opti- 
mization problem (3.3) given by the following partial integro-differential 
equation 



sup (c 7 + ^(t,x,y) 

'0,oo)x[0,ll L ot 



(c,7r)G[0,oo)x[0,l] 

dv 

+ -x-(t,x,y)(TVx(fj,(y) -r(y))+xr(y) - c) 

(3.4) 

1 8 V . . o 2 2/ \ 9v 

+ 2dx 1 ^ t,x,y ' 7T x ° ^ y > ~ Q^^^ x 'y) x y 

+ A / (v(t, x,y + z)- v(t, x, y))v{dz) \ = 0, 

Jz>0 ) 

v(T, x, y) = x 7 . 

As we use a power utility function, it is natural to try to find a solution of 
the form v(t,x,y) = x 7 /(t,y) for some function /. With this choice of value 
function, the optimal strategy (c,tt), which maximizes the left-hand side of 

(3.4) , is given by 

(3.5) c = xf(t,y)- 1 /^\ 

(3.6) vr = argmax{7r(/x(y) - r(y)) - ivr 2 (l - 7)cr 2 (y)}. 

7re[0,l] 

To investigate the formula for the investment strategy more closely, we define 
the three sets 

V 1 = {y>0,^{y)-r(y)<0}, 

(3.7) V 2 = {y> 0,/i(y) - r(y) > 0, (1 - j)a 2 (y) > fx(y) - r(y)}, 
V 3 = {y> 0,/i(y) - r(y) > 0, (1 - j)a 2 {y) < »(y) - r(y)}. 
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The strategy tt is given by 

fO, 

fj>(y) - r (y) 



(3.8) 



TT 



(1- 7)^(2/) ! 



1, 



The following lemma is a counterpart of Lemma 5.1 in [5]. 
Lemma 3.2. Define the function 

1 



Q(y) = max < ir(u(y) — r(y)) 
(3.9) W tt6[o,i]' 



(r(y), 
(V(y)-r(y)) 2 



r(y), 



2(l- 7 )a 2 (y) 

My)- ^0--i)° 2 (y), 



TT 2 {l- 1 )a 2 (y)\+r(y) 

yeVi, 
yeV2, 

yev 3 . 



The function Q is nonnegative, continuous and satisfies the linear growth 
condition 



(3.10) 



0<r(y)<Q(y)<A + By, y > 0, 



for nonnegative A and B. The derivative ofQ is continuous and also satisfies 
a linear growth condition: for nonnegative C and D, we have 



dQ 
dy 



(y) 



<C + Dy, y>0. 



Proof. First, note that the sets T>\ and T>2 have common boundary 

(3.11) dV 12 = {y>0,n(y) = r(y)} 
and that T>i and T> 3 have common boundary 

(3.12) 8V 23 = {y > 0, (1 - 7 )a 2 (y) = fi(y) - r(y)}. 

The sets T>\ and T> 3 do not have a common boundary. 

It is straightforward to show that Q is continuous in T>\, T>2 and T> 3 , as 
well as over the boundaries dT>\2 and dT>23- The linear growth condition 
clearly holds in the sets T>\ and T> 3 . Note that in T>2, the inequality 

(v(y)-r(y)) 2 



(3.13) 



+ r(y)< 1 -(n(y)-r(y)) + r(y) 



2 (My) + r(y)) 



OPTIMAL INVESTMENT AND CONSUMPTION 



9 



holds, from which the linear growth condition of the function Q in the set 
T>2 follows, from (Al). 

We differentiate the function Q and obtain 



dQ 
dy 



(y) 



f dr 

Ty^ 

Ml/) - r{y)){%{y) - %{y)) <M - r(y)) 2 %{y) dr 

(i-7k 2 (y) (i-7k 3 (y) + dy {yh 



dfj, 
< dy 



1 . , da 2 , , 

<y)- 5(1-7)-^), 



Again, it is easy to show that this derivative is continuous in T>i, T>2 and D3 
and over the boundaries dT>\2 and dT>23, and that a linear growth condition 
holds in T>\ and P3. To prove the linear growth condition in the set T>2, note 
that 



dQ 
dy 



(y) 



< 



dfj, 
dy 



(y) 



+ ^(1-7) 



da* 
dy 



(y) 



holds for y € T>2 ■ d 



Remark 3.3. When investigating the unconstrained optimization prob- 
lem, 7r £ 1, the set T>2 must coincide with the whole positive real line 
and one must assume a uniform lower bound of the function a, that is, 
inf ?/> o &{y) > 0. In the Barndorff-Nielsen and Shephard model, this condition 
does not hold unless we introduce reversion to a strictly positive constant 
(i.e., a linear drift term with positive mean reverting level). However, by con- 
sidering a constrained strategy, one can overcome the global lower uniform 
boundedness and work with uniform boundedness only over some subset; 
see [5] and [17] for the structure of the set T>2- Note that for the constrained 
optimization problem, condition (A3) is not necessary. If volatility hits zero, 
one can assume that the set T>2 reduces to an empty set so that the results 
from this paper remain valid. However, in order that all terms in (3.8) and 
(3.9) are well defined, we prefer to retain condition (A3). Moreover, we point 
out that (A3) is very common in stochastic volatility optimization models 
(see [8, 12, 15, 23]) as well as being economically sensible. 

We would like to point out that the linear growth condition (3.10) and the 
relations (2.8)-(2.11) will be frequently applied when proving our results. 
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By substituting (3.5) and (3.6) into (3.4) we arrive at the nonlinear first- 
order partial integro-differential equation for the function /, 

= % (*> f ) - A tt (*> y)y + x l (/(*. y + *) - / (*> j/)M<fe) 

cty Ax) 

(3.14) 

+ if(t, y)Q(y) + (l - y)f(t, y)- 7/(1 " 7) , /(T, y) = 1. 
We will show that there exists a unique classical solution to this equation. 

4. Existence of the solution. We introduce an operator £ acting on func- 
tions / as follows: 



(Cf)(t, y) = E 



3 7 J* QQr*-v(s))da 



(4.1) +(1-7)/ ^^ Q{Yt ' y{ - u))du f{s,Y^{s))-^ l -^ds 



for Q as in Lemma 3.2. By applying (heuristically) the Feynman-Kac for- 
mula to (3.14), we arrive at the following fixed point equation: 

(4.2) (Cf)(t,y) = f(t,y), (t,y) € [0,T] x (0, oo). 

In this section, we prove that equation (4.2) has a unique solution /. In Sec- 
tion 5, we shall show that this solution satisfies the partial integro-differential 
equation (3.14) in the classical sense. 

We start with some observations. Note that it is easy to derive a lower 
bound for the optimal value function V, 

(4.3) V(t,x,y) >xm[e^^ r{Y{s))ds ], 

(t,x,y)€[0,T\ x (0,co) x (0,oo), 

where the left-hand side of (4.3) is the payoff, when the agent does not 
consume and invests everything in the bank account. We conclude that the 
solution to (4.2) should satisfy the inequality 

(4.4) f(t,y)>E[e<^ rins))ds }>l, (t,y)e[0,T\ x (0,oo). 

Moreover, if we apply the operator C to a function / which satisfies (4.4), 
then we obtain a lower bound of the operator, 

(4.5) (jCf)(t,y) >E[e^' r{Y ^ )ds } > 1, (t,y) € [0,T] x (0,co), 

which can be derived by noting that the second term in (4.1) is positive and 
by applying the lower estimate (3.10) of the function Q in the first term. 

We now turn to the more interesting upper bound of the operator C. We 
still assume that condition (4.4) holds, which implies that f(t, y) - T/( 1- T) < 1. 
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By applying the upper estimate (3.10) of the function Q, the estimate (2.11) 
and the representation (2.6), provided that tp^B/X) < oo, we obtain the 
inequality 

(Cf)(t,y) < E L^T-t) +7 B^Y^(s) d s 

+ {l _ l) j\lA(s-t) + ,Br t Y^{u)du ds 
< E e lA(T-i)+CyB/X)y+(yB/X)(L(\T)-L(\i)) 



(4. 6 ) +(1-7) J e fA (s-t)+('YB/\)y+('yB/\)(L(\ S )-L(M)) dfj 

= el A(T-t)+(lB/X)y+X^B/\)(T-t) 

+ (1 - 7 ) [ T e lMs-t)+(lB/\)y+\iP(yB/\)(s-t) dg 



1 - r r \ r A'(T-t)+B'y 
A' J 

where we have introduced the constants A' = jA + Xip(jB/X) > and B' = 
7 5/A>0. 

In the rest of the paper, we assume that the following condition on the 
Levy measure of L, formulated in terms of the characteristic exponent in 
(2.6) holds: 

(B) < 00 for w = 2(1 + l)(B' V B' a ) + e and some e > 0, 

where B' a = ^B^/X > is defined analogously to B' and B a is defined in 
(Al). The reason for this assumption becomes clear in the course of our cal- 
culations. It is needed in Section 6 in order to verify the optimality. Note that 
the lemmas in this section and Section 5 hold under integrability conditions 
of lower orders. 

Let us investigate the operator C in a more rigorous way. Denote by 
C e ([0,T] x (0,oo)) the space of continuous functions / on [0,T] x (0, 00) 
satisfying 

1 < f(t, y)<(l + ~~xr~} e A ' (T -' )+B \ (t, y) G [0, T] x (0, 00). 



A' 

We define a metric an C e ([0, T] x (0, 00)) by 

(4.7) d(<p,0= sup \e- a ^- B 'y^(t,y)-at,y))\, 

(t,2/)e[0,T]x(0,oo) 

for some a > A' to be specified later. The space (C e ([0, T] x (0, 00)), d) is a 
complete metric space. Below, we state two lemmas dealing with the prop- 
erties of the operator C. 
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Lemma 4.1. The operator C defines a mapping from C e ([0,T] x (0, oo)) 
into itself. 

Proof. Based on our previous results (4.5) and (4.6), we can conclude 
that the lower and upper bounds are preserved. It remains to prove the 
continuity of the mapping (t,y) t— > (Cf)(t,y). Due to the time homogeneity 
of Y, the operator C can be represented as 



,1 fl * QO^'H'))^ 



+ (1 - 7 ) F* e^o Q(y°^))du f[s + ^ F o )2/(s)) -7/(i-7) ds 



(Cf)(t,y)=K 



The above representation simplifies proving continuity in the time variable. 
Note that by the growth condition (3.10) and relation (2.11), 

g 7 /„' " + ^ y0,«/ a ^-7/(l-7) < g 7 /„' Q(y '»W)efa 

< e 7AT+B'3/+S'L(AT) 

holds P-a.s. and the cadlag mapping (y, u) i— > Y 0,?/ (u) is bounded a.s on com- 
pact sets. In order to prove continuity in the time variable, one can directly 
apply Lebesgue's dominated convergence theorem and take the limit under 
the integral. To prove continuity of the mapping y \— > (Cf)(t,y) at a fixed 
point yo > 0, define a compact set J7 around yo an d take a sequence of points 
y n EU such that y n — > yo as — > oo. In this setting, we can find a uniform 
bound for all y n and can apply Lebesgue's dominated convergence the- 
orem. The continuity of (t, y) i— > (Cf)(t,y) now follows from the continuity 
of / and Q and the continuity of the mapping y h- » Y '^. □ 

Lemma 4.2. Tfte mapping £:C e ([0,T] x (0,oo)) ->C e ([0,T] x (0,oo)) zs 
a contraction with respect to the metric (4-7) for a > A' + 7. 

PROOF. Take two functions </?,£€ C e ([0,T] x (0, 00)). Again, we invoke 
(2.10) and (3.10). The following inequalities then hold for all (t,y) € [0,T] x 
(0,oo): 

d(£<p,£Z) = |e- Q ( T -')^((£^)(t,y) - (^)(t,y))| 

r .7/;Q(^' ! 'H)^| v9 ( S) yt, ? /( s ))-7/(l-7) 



X E 



T 

er ' 



■ c(s,y i,y (s))" 7/{1 " 7) Ms 
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< le -a{T-t)-B'v E 



?f*' Q <r t,v w* t \< P (8,Y t >y(8)) 

-£( S ,y*"( S ))|cfe 



-a(T-t)-B'y 



3 7 J t s Q(Y t 'y(u))du+ a {T-s)+B'Y t 'y(s) ^ 



(■T 



x E 



■yA(s-t)+B'{y+L(\s)-L{Xt)~Y t -y{s))+a(T-s)+B'Y t 'y(s) 



-a(s-t)+yA(s-t)+\il)(B')(s-t) dg 



ds 



< 



1 



a — A' 

where the mean value theorem has been applied in line 3. We conclude that 

dOC</>,££)<Cd(</?,£), C<1, 
which proves that the operator C defines a contraction mapping. □ 

The main result of this section is the following proposition, which is a 
consequence of Banach's fixed point theorem. 

Proposition 4.3. The equation 
(4.8) (Cf)(t,y) = f(t,y) 

has a unique solution /€C e ([0,T] x (0, oo)). 

5. Differentiability of the solution. In this section, we establish the dif- 
ferentiability of the function /. In order to apply a classical verification 
theorem, we have to prove that / is continuously differentiable in the time 
and in the space variable. 

We assume that B > 0. In the case when the function Q is uniformly 
bounded in y, that is, B = 0, an arbitrary, strictly positive (small) constant 
B > can be chosen so that the proofs from this section remain true. We 
remark that our arguments can be modified in order to handle the special 
case of B = and to derive sharper bounds. We would like to point out that 
the main theorem of our paper, Theorem 6.1, holds true even for B = 0. 

Recall that the ODE 



(5.1) 



^t(t) + ( 1 -\)<j>(t) + \a = 0, (f>(T) = a, 
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has the unique, smooth and strictly positive solution in the class C 1 ([0,T]) 
given by 

0(i) = a + 7 / 0(s)e- A(s_t) ds, 



with constant a > 0. 

The idea for establishing differentiability in the space variable is to con- 
struct a sequence of functions (/ n )neN which converge to / and which share 
some desirable properties. 

Lemma 5.1. Define A" = -yA + \^(B") > and B" = B'(l + J) > 0, 
with B' = -yB/X > and A,B as in (3.10). 

Choose a function ft eC e ([0,T] x (0, oo)) !~l C 0,1 ([0, T] x (0,oo)) such that 



(5.2) 



df 



dy 



-fry) <cp(t)e A "^ +B "y, (t,y) € [0,T] x (0,oo), 



where 4> solves (5.1) with a = j(l + ^7r)(^r V C7) > and C,D as m 
Lemma 3.2. 

Now, construct now the sequence (f n )neN recursively as f n +i = £fn with 
£ defined as in (4-1)- 
Then, for all n G N, 

/ n GC e ([0,T] x (0,c 5 o))nC°' 1 ([0,T] x (0,oo)) 



and 
(5.3) 



< *(() e - 4 "< r -" +B "», (t,») e[o,r]x (0,00). 



Proof. Recall from (4.1) that 
f2(t,y)=E 



+ (i- 7 ) /' T e ^/Q( Yt ' ! '( u )) du / 1 ( s ,y**( s ))-T/(i-T) ds 



First, we prove that the mapping (t,y) 1— > ^-(t,y) is continuous. 



We expect that the derivative equals 



9£ 
<9y 



(t,y)=E 



7e 7j? / ^(Y*.V( a )) e -A(-*) da 



+ 



(1-7) / T 7 e^; ( 3( yt ' !/ ^/ 1 ( s ,y( s ))^/( 1 ^) 



(5.4) 



dQ 



(y t ' J/ (n))e" A(,1 - i) (iuds 
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' 7e 7/;Q(y t -w)^ /l(S)yt ,y (s)) -i/(i-7) 



oy 

This will follow from Lebesgue's dominated convergence theorem, provided 
that it can be applied. Below, we establish three estimates which allow us 
to interchange differentiation and integration. We point out that the inter- 
change is justified if we can bound the derivative by an integrable function. 
The estimates are also used later to establish (5.3). 

We recall that in order to find a uniform bound, one can take a limit 
Vn —> 2/0) as n — > oo, over a sequence of points y n € U, where U is a compact 
set around a fixed point yo > 0. 

Note that by invoking the simple inequality a + by < (- V a)e bey for all 
e > 0, together with (2.9), we obtain that 



(5.5) 



< (C + DY t >y{u))e- x{ - u -^ 
4D 



< 



< 



\lB' 



V C ) (JaB' /4){y+L{\s)-L{\t)) e ~X(u-t) 



holds P-a.s., for <t <u< s <T. Based on (5.5), we derive that 

— ( e ~< J t s Q(Y t ' y (u))d^ 
dy 



le ,J:Q(Y^(u) d ud_ ^ Q{ yt,y {u)) du 



(5.6) 



1 e 



7 //Q(y*'»)(u)du 



d_ 

dy 



{Q{Y^{u)))du 



< 



x e -yA(s-t)+B'(y+L(Xs)-L(Xt)) e (-yB'/4:)(y+L(Xs)-L(Xt)) 
< I f^R V Cl \ e lA(s-t)+B"{y+L(Xs)~L{X£)) 



-X(u-t) du 



\\B> 

holds P-a.s., for < t < s < T. We would like to point out that we are allowed 
to interchange integration and differentiation in the first line of (5.6) since 
the bound (5.5) is integrable P-a.s. on [t,s]. 
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Based on (5.2) and (5.6), we obtain the third estimate 
d 



dy 



((l- 7 ) e 7/;Q(^(«))^ /l(S) y (s)) -7/(l-7) ; 



(1 _ 7) ^ (e 7/;^(^ )/l(S)y(s)) -7/(i-7) 

< ^~ 7 ^ ( — V C 7 ^) e 7A(*-t)+B"(y+i(A S )-L(At)) 



(5.7) +7e' 



-A(s-t) 



< 



(1-7) /4D 



J V C 7 ^) e 7^(s-i)+B"(y+i(A S )-L(At)) 



A V# 

+ 7e 7^(^-t)+fi''fe+L(A S )-L(At)-y*^( S ))^^^A''(T- S )+B''y t .n S ) e -A( S -t) 



1-7 /4D 



37 A(s-t)+B"(y+L(As)-L(At)) 



'-a.s. 



+ 7 ^( s ) e -A(«-i) e A"(T-s)+ 7j 4( S -t)+B"j/+B"(L(A S )-L(At)) ^ 

As the derived bound (5.7) is a cadlag mapping, it is a.s. integrable, and we 
have that 

8 i-t 



dy 



((l_ 7)e 7/;Q(^H)^ /l(S)y(s)) -7/(l-7) )ds 



T d_ 

dy 



( (i_ 7)e 7/;Q(^w)^ /l(S)y(s)) -7/(i-7) )dS5 



-a.s. 



Finally, taking the derivative under the expectation is also justified since, 
by condition (B), we have ip{B") < oo. Consequently, we have shown that 
the derivative (5.4) holds. 

The continuity of the mapping (t, y) i— ► ^-{t, y) again follows from Lebesgue's 
dominated convergence theorem, by applying the estimates (5.6) and (5.7), 
and from the continuity of the functions /i and Q, as well as their derivatives 
[cf. the proof of continuity in Lemma (4.1)]. 

We still have to prove that the bound (5.3) holds. By combining (5.6) and 
(5.7), we can estimate for n = 2: 



1 /AD 
\\~B 



j V OV) e lMT-t)+B"(y+L(\T)-L(\t)) 
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+ 



1-7/4D 



A V B 



7 VC7 7 



T 



^A(s-t)+B"(y+L(Xs)-L(Xt)) dg 



+ 7 



T 



^ e -X(s-t) e A"(T-s)+~fA(s-t)+B"(y+L(Xs)-L(\t)) dg 



Lf^E. yCi yA"(T-t) + B»y 



X\B 

1 -7 {AD 



+ 



A V B 



7 vc 7 



r 



+ e A"(T-t)+B»y 7 J- r (f){s)e -X { s~t) ds 



< 



1 /4L> 



7 VC7 7 )e A "( T - i ) +B "^ 



A \B 

1-1 (AD 



+ 



\A" V B 



,A"{T-t)+B"y 



+ e A"{T~t)+B" yi 



[s e 



-X(s—t) ds 



A"(T-t)+B"y 



a + 7 



5 )e- A ( s -*)d S )=^(t)e A "( T - t )+ s "^ I 



where we have invoked the solution of the ODE (5.1) with the appropriate 
constant. Repeating the calculations recursively concludes the proof. □ 

From the properties of the constructed sequence (f n )neN, we can deduce 
an important property of the function /. 



Proposition 5.2. The function f belongs to the class C e ([0, T] x (0, oo)) n 
C 0,1 ([0,T] x (0, oo)). Moreover, its derivative satisfies 



(5.8) 



A"{T-t)+B"y 



(t,y)€[0,T\ x (0,oo), 



for A" , B" and <fi as in Lemma 5.1. 



Proof. The result follows if we show that the sequence (^ L ) n eN ; con- 
structed in Lemma 5.1, converges uniformly, at least on compact subsets of 
[0,T] x(0,oo). 

Choose n>m and p> a V A". Using the definition of the derivative (5.4), 
we have 



-p(T-t)-2B"y 



dfn+1 (t,y)-^(t,y) 



9y 



9y 
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(i_ 7 )y%7/;Q(^n«))^ 

-/ m ( s ,y^( s ))-7/(i-7)| 

A ^-*)(u) duds 



g<r*M) 



-p(T-t)-2B"y 



+ E 



(5.9) 



y T 7e 7rQ(^-w)^| /n(s ,y*,i/ (S ))-i/(i-7) 

-/ m (^(*)r 1/(1 ~ 7) l 

x ^( Si ^( s )) e -^-') da" 



-p(T-t)-2B"y 



+ E 



r 



7 e 



7 //0(y*-»(u))d« 

^( s ,y^ (s ))-^( S) y^( s )) 



X /„(S, y^( S ))- 1 /(l-7) e -A( S -t) ds 



-p(T-t)-2B"y 



=:Mi + M 2 + M 3 . 

We first derive an upper bound for Mi. Again, let d(-, •) denote the metric 
defined in (4.7). Applying the estimate (5.5) and the mean value theorem, 
we find 

rT 



Mi < 7^E 



( T/A(s-t)+B'(y+L(Xs)-L(Xt)-Y t -y(s)) e p(T-s)+B'Y t 'y(s) 



x e 



t(T-s)-B'Y t 'y(s) 



(5.10) 



x e 



x|/ n ( s ,y^( s ))-/ m ( s ,y^( s ))| 

X ^ vC ) e hB> my+ L(Xs)-L(Xt)) l~ e -Hu-t) duds 
-p{T-t)~2B"y 



< 



1 f4Dj 
A 



( ^vc 7 2 jd(/„,/ m) y e (7A + A* ( B«)-p)(-o ds 



<^ld(/n,/m)- 
Similarly, we have 



r 



Mo < 



1-7 



-E 



e7j 4( s -i)+B"(y+L(A S )-L(At)-y''!'(s)) e p(T- S )+B'y i ^( S ) 
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x e -^-)- s '^W| /n ( s ,y^( s ))_ /m(s ,y^( s ))| 



x e 



1 



-p(T-t)-2B"y 



(5.11) <-^d(/ n ,/ m )E 

1 - 7 



D jA(s-t)+B"(y+L(Xs)-L(\t)-Y t 'y(s)) 



x e 



p(t-s)+S'(2/+L(As)-L(At)) 



x<p(s)e A "^ +B " Yt ' v ^e- x ^e- 2B "yds 



< ^(U, f m )e A " T sup {#*)}/ 



i€[0,T] 



<K 2 d(f n ,f m ), 



where we have used the bound (5.3) for the sequence of derivatives (^r)neN- 



Finally, we obtain a bound for M3: 

i-T 



<E 



(5.12) 



7e 



x e 



7A(s-t)+2S"(3/+L(As)-L(At)-y''»(s)) 



-p(T-s)-2B"y t -«(s) 



x eP{ T-s) + 2B"Y^(s) e -X(s-t) dg 



x e 



-p(T-t)-2B"y 



< sup 

(t,3/)e[0,T]x(0,oo) 



-p(T-t)-2B"y ( 9fn 



7/ e (7A+AV>(2B'')-p-A)(5-t) ^ 



< K3 sup 

(t,j/)G[0,T]x(0,oo) 



-p(T-t)-2B"y 



where p must be chosen such that K$ = j(p — jA — Xip(2B") + A) -1 < 1. 

Note that by the contraction property of the operator C proved in Lemma 
4.2, we have 

m— 1 



d(/n,/m) < 



(5.13) 



1 



<2 



a- A' 

7 

a- A' 



d(/ n -m+l, fl 
"»"!/ I-7 
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By combining (5.10)-(5.13), we get 



-p(T-t)-2B"yf dfn+1 



<2(K 1 +K 2 ) 
+ K 3 



7 



a- A' 



m—l 



1 + 



1-7 
A' 



sup 

(i,3/)G[0,T]x(0,oo) 



e -PCT-*)-^ w( ^ (tty) _^ (tjy) 



< 2(^+^2) 



7 



a- A' 



m—l 



1 + 



l-^l-i^" 1 



A' 



1-K S 



+ Kip' 1 sup 

(t,j/)e[0,T]x(0,oo) 



-P(T-t)-2B''y^f^ M _^ y) 



<2{K 1 +K 2 ) 



7 



m—l 



1 + 



1 - 7 \ 1 - if. 



A' 



m—l 



■ + 21T™- 1 sup {0(i)}, 
1 _ A 3 te[o,T] 



from which we conclude that the sequence (-^-(i, y))neN converges uniformly 
on compact sets. □ 



We now turn to the question of differentiability in the time variable. We 
first show that the function f(t,y) belongs, for every fixed t 6 [0,T], to the 
domain of the infinitesimal generator of the process Y; see Chapter 1.3 in 
[19]. 

As the mapping y i— > f(t,y) is continuously differentiable on (0, oo), we 
can apply Ito's formula and show that the limit relation 



(5.14) 



lim 



E[f(t,Y°<y(s))]-f(t,y) 



-^(t, y)Xy + \ (f(t, y + z)- f(t, y))v{dz) 



2>0 



holds, provided that, for s > 0, 



(5.15) E 



Jz>0 



(f(t,Y°>y(u-) + z)- f(t,Y°>y(u-)))N(du x dz) 



0. 



where N(du x dz) := N(du x dz) — v{dz) du is the compensated Poisson 
random measure from (2.5). It is well known (see, e.g., Theorem 4.2.3 in [1]) 
that condition (5.15) is equivalent to 



(5.16) E 



Jz>0 



\f(t, Y°>y(u-) +z)- f(t, Y^(u-))\ 2 v{dz) du 



< oo. 
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E 



(5.17) 



Jz>0 

<E 
< E 



\f(t, Y°>y(u-) + z)- f(t, Y^{u-))\ 2 v{dz) du 
^ {t yA" { T-t) + 2B"(y^ { u) +z)z 2 v{dz)du 

2 / . \ 2A" (T-t)+2B"y+2B"L(XT) 2B"z 2 



Jz>0 



Jz>0 



P{t)e 



e lu z z l v{dz) du 



< Ke 2B " / z 2 v(dz) + K / e 2B " z z 2 u{dz) 

J0<z<l Jz>\ 

for some positive constant K, which is finite since ip(2B") < oo. The first 
term in (5.17) is clearly finite. We show that the second term is also finite. 

By applying the inequality z < -^e 4 2 and assumption (B), we find that 



{2B" +1 B'/2)z dz 



jB'J J z >! 



,2B' (1+7/2)2 dz 



< OO. 



Before stating the next lemma, we would like to remark that the mapping 
(t,y) 1— > f z>0 (f(t,y + z) — f(t,y))u(dz) is continuous on [0,T] x (0,oo). This 
follows from the inequality 

\f(t, y + z) - f(t, y)\ < me A " {T ~ t)+B " {y+z) z 

and Lebesgue's dominated convergence theorem. 

Proposition 5.3. The function f satisfies the partial integro- differential 
equation 



(5.18) 



at ay J z> o 

+ f(t,yhQ(y) + (1 - 7 )(/(t,y))~ 7/(1 ~ 7) , f(T,y) = 1, 



in the classical sense. In particular, the mapping (t,y) 1— > ^£(t,y) is contin- 
uous on [0,T) x (0,oo). 



Proof. The idea of the proof is similar to that of the proof of Propo- 
sition 5.5 in [5]. We will calculate the limit in (5.14) explicitly by using the 
representation of /. 
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Consider a fixed t € [0, T). Note that by the time homogeneity of Y, the 
equivaient representation holds: 



f(t,y) :=E 



T-t+s 



Q(K(u;))dio 
T-t+s 



+ (1-7) 



7 /" Q(r(™))du) 



e J « 



x (/(u + t-5,y(«))) 



-7/(1-7) 



du | V(s) = y 



for s > 0. Let o"((y 0,y (s)) denote the cr-algebra generated by the random 
variable Y°' y (s) as defined in (2.7). We have that 



f(t,Y°>v(s)) 
= E 



e 7 I 
+ (1-7) 



Q{Y{w))dw 
T-t+s 



7 £Q(Y(w))dw 



e J s 



x (/(« + 1 - a, y(n)))- 7/(1 " 7) | <r((y°'f (s)) 



holds P-a.s. Applying the law of iterated expectations, we obtain 

E[/(t,y°-»( s ))] 



E 



e 7 J S T t+s Q(y(»))rf«- 



+ (1-7) 



e 7 /"0(y(iB))dui 

x (/(« + 1 - s , y(u)))" 7/{1 " 7) du I y(0) = y 



Now, consider the difference K[f(t,Y 0,y (s))] — f(t,y) for some s > in the 
neighborhood of 0. By simple algebraic manipulations, we find 



- s (E[f(t,Y°'y(s))]-f(t,y)) 



E 



(1-7) 



T-t+s 



( / (u + t _ S) y0,, (u))) -7/(i-7) e 7/;Q(^M)^ 
xi( e -^/o*«t^"(«))*»_i )du 



-E 



(l- 7 ) / (/(n + t- S ,y°^(u)))~ 7/(1 ~ 7) e 7 /o" ( 3( y0 ' !/ ^dn 
o 



OPTIMAL INVESTMENT AND CONSUMPTION 



23 



+ E 



{i Jo tf "<3(r° iy M)<ft*I( e -7 f*Q(Y°>y(w))dw _ y 



+ l(f(t-8,y)-f(t,y)) 
s 

=: Mi(s) + M 2 {s) + M 3 (s) + M A (s). 
Note that 

g 7 f*Q(Y°'«(w))dw < ^yAT+B'y+B' L(XT) 

— (1 - e-T/o'Q^^W)^) < sup{7 Q(yO,, (s))} 

<7^T + 7 J Bsup{y '' i/ (s)} 

s>0 

< 7^T + 7£y + jBL(XT) 

<(\y 1 AT)e B 'y +B ' L ^ 

hold P-a.s. for < s < T. The above estimates ensure that we can apply 
Lebesgue's dominated convergence theorem to obtain the following limits: 

limMi( S ) = -7Q(y)E[(l- 7 ) F'* #&QV*' V W)** 
s ^o L Jo 

x (f(u + t,Y°'y(u)))' 7/{1 - 7) du 

hmA/ 2 ( S ) = -(l- 7 )(/>,y))" 7/(1 - 7) , 

limM 3 (s) = - 7 Q(y)E[e 7 /o T " i( 3( y0 ' ! '(-))^]. 

Moreover, lim s _>o(jMi('S) + M^{s)) = —jQ(y)f(t,y) holds and by combining 
these calculations with (5.14), we arrive at 

hmM 4 (,) = -lim /( ^ ) - /(t - g ' y) 



?f(t,y)\y + [ (f(t,y + z)-f(t,y))v(dz) 
oy Jz>o 



dy 

+ f(t,y)iQ(y) + (i—f)(f(t,y)r i/{1 " r) - 

We conclude that the derivative 4j£ exists and that / satisfies the partial 

integro-differential equation (5.18). Moreover, the mapping (t, y) i— > ^{t,y) 
is continuous on [0, T) x (0, oo) by the continuity of all terms on the right- 
hand side of (5.18). □ 
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We can also conclude that the function / is the only classical solution of 
the partial integro-differential equation (3.14) as, for any such solution, the 
Feynman-Kac representation must hold; see [4] for a similar argument. 

6. Optimality of the solution. We shall conclude with the following the- 
orem, which states that our solution is indeed optimal. 

Theorem 6.1. Assume that the conditions (A1)-(A3) and (B) hold. 
Define the investment strategy 

(6.1) 7r(t) = arg max {7r0x(y(t-))-r(y(t-)))-|7r 2 (l- 7 )a 2 (y(t-))} 

7re[0,IJ 

and the consumption rate 

(6.2) g(t) = x a '* (*)(/(*, r(H)r 1/(1 ~ 7) , 

where the function f is the unique solution of the fixed point equation (4-8) 
in the space C^flO.T) x (0,oo)) PlC e ([0,T] x (0,oo)) given by 



3 *tf t Q(X*'Hs))d8 

+ (1 - 7 ) jf ft*' Q0f*Hu))du f ^ Y t,y {s)) -,/(i-,) ds 

and X c ' w is the wealth process of the agent under (c,7r), defined as 

X ^(t) =xe / t ^( s )^( y ( s -»- r ( y ( s -))) +r ( y ( s -))-^ s ' y ( s -))^ 1/(1 ~ 7) ) ds 
(6.3) " 

x e -l/2 j*(ir(s))2^(Y(s-))ds+f*n(sMY(s-))dW(s)_ 

The pair (c, n) is then the optimal strategy for the investment and consump- 
tion problem (3.2). 

The proof of the above theorem is based on a verification theorem stating 
the conditions which a candidate value function should satisfy in order to 
coincide with the optimal value function. Theorem 3.1 of [19] is an appropri- 
ate verification theorem for jump-diffusion processes. We can prove that our 
candidate solution satisfies all of its conditions. In particular, we would like 
to point out that assumption (B) is needed for proving uniform integrability 
of the value function. We want to emphasize that we have been able to prove 
the optimality of the strategy under the weaker integrability assumption (B) 
on the Levy measure of L than is required in [5] . 

Finally, we want to mention that in [11], a different verification theorem 
for jump-diffusion processes has been proven, one which requires substan- 
tially weaker conditions than those in Theorem 3.1 of [19]. However, this 
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result requires the strategy to be caglad (left continuous with right lim- 
its), in contrast to the weaker progressively measurable condition as in the 
present paper. Details can be obtained from the authors on request. 

Remark 6.2. Assume that we state our problem for a deterministic 
function Y. The candidate for the value function is again V(t, x) = x 7 /(t) 
resulting in the ODE 

f (t) + f{t)lQ(Y{t)) + (1 - 7)(/(i))" 7/(1 " 7) = 0, f(T) = 1, 
which has the solution given by the fixed point equation 

jit) = cfttfQinsVds + (1 _ 7) r T e ~/j;Q(Y(^ f - {sr y/(i^) ds , 



0<t<T. 

If Y is a stochastic process, then the function /, defined as above, involves 
a random path of the volatility process, so it depends on u € ft. It is tempt- 
ing to believe that the function E[f(t,u) \ Y(t) = y] solves the optimization 
problem (3.2). 

However, taking the operator C as defined in (4.1), we calculate 



= E t ' y 



<E f ' y 



+ (1 - 7 )^ T e 7 r Q{y( ^ ))du (E[/( S ,o;)|^ s ])-^( 1 -^ ds 

e 7 J t T Q(Y(s,u))ds 
T 



+ (l- 7 ) J e^t Q{Y ^ ))du f(s,u)-~<^ l -~<Us 
= E t >v[f(t,u)}, 

where the equality holds if and only if (Y(t,u>))o<t<T is independent of to, 
hence deterministic. 

We conclude that the function E* ,y [/(t,u;)] does not satisfy the fixed point 
equation (4.8) and, as a result, it is not the solution to our optimization 
problem. The optimal value function and the optimal investment and con- 
sumption strategy are, different, as one might have expected. 



Example 6.3. We consider a financial model of Barndorff-Nielsen and 
Shephard type. More precisely, we choose the time horizon T = 1 and 7 = 
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Fig. 1. The optimal consumption rate as a function of time to maturity and the volatility 
level; see also the text in Example 6.3. 

0.75 for the exponent of the power utility function. Furthermore, in (2.1), 
we choose A = 1/6, in (2.2), we take r(y) = and in (2.3), we take fx(y) = 

0. 1 + 0.5y and <J 2 (y) = y. Let the subordinator L be a compound Poisson 
process with jumps of intensity 0.5 and exponentially distributed jump sizes 
with expectation 1/15. We set the initial volatility level at Y(0) = 0.2, which 
equals the expected long-term volatility. 

We have solved the nonlinear partial integro-differential equation (3.14) 
numerically by applying an explicit finite difference method. As we are deal- 
ing with a first-order integro-differential equation and the Levy measure is 
finite, the explicit scheme is more efficient than the implicit scheme; see [9], 
Chapter 12.4, for details. We point out that the finite difference method has 
been applied to the transformed equation to which the solution is f(t, y)e~ Ky . 
The exponential scaling has been applied in order to set a sensible bound- 
ary condition in the bounded domain. The parameter k should be chosen 
sufficiently large so that lim^oo f(t,y)e~ Ky = holds. 

Based on (3.8), we can state that the optimal investment strategy is 7r(t) = 

1, whereas the optimal consumption rate is given in Figure 1. 

The first observation, which is common in optimal investment and con- 
sumption models, is that the optimal consumption rate is an increasing 
function of time. In the model considered, it is interesting to note that the 
optimal consumption rate is a decreasing function of volatility level. The re- 
sult agrees with our intuition: the higher level of volatility leads to a higher 
variability, which is, however, compensated for generously by an increase 
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Fig. 2. The optimal consumption rate in the stochastic volatility model (upper curve) 
and in the model with constant volatility (lower curve); see also the text in Example 6.3. 

in the appreciation rate of the risky asset. This explains why the investor 
should consume less and invest more. 

It has already been stated, in [5], that stochastic volatility modeling can 
change the investment strategy significantly. We have simulated one path 
resulting in a volatility that jumps at t = 0.05 by 0.12 and at t = 0.65 by 
0.07. The optimal consumption pattern is significantly different, when com- 
pared with the constant volatility model Y = 0.2; see Figure 2. Under the 
stochastic volatility model, it is optimal to consume much higher propor- 
tions of the wealth as the unexpected jump in the volatility increases the 
variability of the return and may cause a severe decrease in the portfolio 
value. Note the discontinuity in the consumption strategy at t = 0.05 in the 
upper curve in Figure 2, which is caused by the jump in the volatility. The 
second discontinuity at t = 0.65 is much less visible. 

We conclude this paper with the solution for the optimization problem in 
the case of a logarithmic utility. Logarithmic utility is investigated in depths 
in [14], where a general financial market is considered, consisting of stocks 
whose prices are driven by semimartingales. The solution is stated in terms 
of the semimartingales, characteristics, which are not straightforward to find 
in our model. 

We make an ansatz with a value function of the form 
v(t,x,y) = g(t) log x + h(t,y). 
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This yields the following equations: 

(6.4) ° = ^ + 1 ' 9 ^ = 1 
and 

dh . . dh , 

(6.5) +\[ (h(t,y + z)-h(t,y))u(dz) 

Jz>0 



where 



+ g(t)Q°(y) - logfl(t) -i = o, h(T, y) = o, 



Q°(y) = max {vr(/i(y) - r(y)) - ±vrV(y)} + r(y) 
7re[0,lJ 



is the analogue of (3.9). 

The following theorem can be proven. 

Theorem 6.4. Assume that conditions (A1)-(A3) hold and that the 
Levy measure v of L satisfies the following condition: 

(C) J Z>1 z 1+£ v(dz) < oo, equivalently E[L(l) 1+e ] < oo, for some e > 0. 

Define the investment strategy 

7r(t) = argmax{7r(/i(Y(t-)) - r(Y(t-))) - \n 2 a 2 {Y(t-))} 
7re[0,l] 

and the consumption rate 

X £ ^(t) 



m 



l + T-t' 

where X C,7T is the wealth process of the agent under (c,fr), defined as 

X & '*U) = x Jo(Hs)MY(s-))-r{Y(s-)))+r(Y( S -))-l/(l+T-s))ds 
x e -l/2 J*(Hs)) 2 a 2 (Y{s-))ds+J*7t(s)a(Y(s-))dW(s)_ 

The pair (c, fr) is then the optimal strategy for the investment and consump- 
tion problem under a logarithmic utility function. 

As we are facing the linear equation (6.5), existence and smoothness of a 
solution can be easily proven by combining the results from this paper with 
those and from [5]. It is well known (see [18]) that the optimal consumption 
rate in the case of a logarithmic utility does not depend on the financial 
coefficients. 
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7. Conclusions. In this paper, we have solved an investment and con- 
sumption problem for an agent who invests in a Black-Scholes market with 
stochastic coefficients driven by a non-Gaussian Ornstein-Uhlenbeck pro- 
cess. We have proven that the candidate value function is the classical so- 
lution of the corresponding Hamilton-Jacobi-Bellman equation. In partic- 
ular, we have provided a classical solution to a nonlinear first-order partial 
integro-differential equation. 

The optimal investment strategy has been explicitly calculated, while the 
optimal consumption rate depends on the function which solves the partial 
integro-differential equation. The conclusion from the simulation study is 
that under stochastic volatility, the optimal consumption strategy is signif- 
icantly different compared to a constant volatility model. 

In [5], a multivariate Ornstein-Uhlenbeck process driven by independent 
subordinators was considered, while in [17], a financial market consisting of 
n stocks was investigated. We would like to point out that our results can 
be extended to both settings. 
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